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Abstract 

Numerical  modeling  of  the  evaporation  process  in  sprays  under  diesel  conditions  is  key  for  the  development  of 
efficient  injection  strategies  and  to  increase  combustion  efficiency.  In  this  study,  three-dimensional  numerical 
simulations  of  single  and  two-hole  injector  nozzles  under  diesel  conditions  are  conducted  to  study  the  spray 
behavior  and  the  effect  of  multi-hole  nozzles  on  heavy  fueled  spray  parameters  and  mixing.  The  configuration 
corresponds  to  a  high-pressure  JP-8  spray  injected  into  a  high  temperature  pressure  vessel  (HTPV)  flow-through 
combustion  chamber  simulating  realistic  conditions  found  in  diesel  engines.  A  Lagrangian  particle  tracking  method 
coupled  with  a  classical  blob  injection  wave-based  model  is  adopted  through  the  use  of  CONVERGE  solver  to  treat 
the  spray  process.  An  Adaptive  Mesh  Refinement  (AMR)  and  fixed  embedding  technique  is  employed  to  capture  the 
gas-liquid  interface  with  high  fidelity  while  keeping  the  cell-count  reasonable.  Two  JP-8  liquid  fuel  surrogates, 
Surrogate-C  (of  60%  n-dodecane,  20%  methylcyclohexane,  20%  o-xylene),  and  the  Modified  Aachen  surrogate 
(80%  n-dodecane  and  20%  trimethylbenze)  were  interrogated  resulting  in  good  agreement  with  high  fidelity  spray 
measurements.  Spray  simulation  results  are  compared  to  our  in-house  experimental  data  for  JP-8  with  single  axial 
hole  and  two-hole  adjacent  (60°)  nozzle  configurations.  Standardized  Engine  Combustion  Network  (ECN)  Spray  A 
ambient  conditions,  consisting  of  900  K  and  60  bar,  are  selected  with  a  rail  pressure  of  1000  bar  with  nominal 
nozzle  diameter  dj  —  147  jim.  Utilizing  k  —  e  Reynolds-Average  Navier  Stokes,  and  dynamic  structure  Large  Eddy 
Simulation  methodologies  both  configurations  result  in  a  20  mm  mean  liquid  length  for  both  single  and  two  hole 
cases.  This  is  in  good  agreement  with  experiments  with  the  single,  while  for  the  double  hole  this  is  an  over¬ 
prediction  of  3-5  mm.  The  reported  differences  are  partially  attributed  to  the  internal  nozzle  flow  dynamics,  reported 
nozzle  eccentricities,  the  influence  of  these  effects  are  discussed. 
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Introduction 

In  many  combustion  systems,  fuel  atomization  and 
the  spray  breakup  processes  play  an  important  role  in 
determining  combustion  characteristics  and  emissions 
formation.  Due  to  the  ever  rising  need  for  greater  fuel 
efficiency  and  lower  emissions,  the  development  of  a 
fundamental  understanding  of  its  process  is  essential 
and  remains  a  challenging  task.  In  U.S.  military 
unmanned  aerial  system  (UAS)  and  unmanned  ground 
system  (UGS)  applications,  the  characterization  of 
engines  operating  with  JP-8  aviation  fuel  is  vital  to 
understand  engine  performance,  combustion  phasing, 
and  emissions  when  JP-8  is  fully  substituted  for  diesel 
fuel  [1].  This  emerges  from  the  Department  of  Defense 
(DOD)’s  “Single  Fuel  Forward”  policy  requiring 
military  vehicles,  and  equipment  to  be  capable  of  using 
a  single  battlefield  fuel  [1].  Hence,  characterization  of 
JP-8  spray  and  combustion  will  remain  of  critical 
priority  in  enabling  the  Army  mission  while  facilitating 
fuel  logistics  operations. 

In  recent  years,  the  engine  modeling  community 
has  relied  on  the  Engine  Combustion  Network  (ECN) 
as  the  premier  research  consortium  to  validate 
physical/chemical  models  and  demonstrate  several 
numerical  features  of  CFD  solvers.  ECN  group  includes 
multi-institutional  facilities  leveraging  expertise  and 
characterizing  diesel  and  gasoline  sprays  through  use  of 
optically  accessible,  constant-volume  combustion 
vessels  under  simulated,  quiescent  diesel  engine 
conditions  [2-3].  Much  of  the  modeling  work  is  aimed 
at  demonstrating  features  important  for  multiphase  3D- 
CFD  solvers  such  as:  the  ability  to  accurately  describe 
Lagrangian  particles  in  Eularian  gas  phase  models, 
assess  the  performance  of  leading  turbulence  models, 
model  combustion  through  adequate  surrogate 
selection;  and  numerically  to  obtain  grid  convergent 
simulations  with  low-to-high  fidelity  models  of  global 
spray  and  combustion  properties  [4],  [5]. 

Using  a  RANS  technique,  Som  et  al  [5]  conducted 
several  validation  studies  and  presented  an  assessment 
of  leading  Lagrangian/Eularian  atomization  models  in 
diesel  non-reacting  sprays.  Some  discrepancies  in  the 
atomization  models  were  reported  for  non-evaporating 
spray  conditions  in  terms  of  the  dispersion  behavior 
while  this  effect  was  not  as  pronounced  in  evaporating 
sprays  where  good  agreements  were  obtained.  It  also 
was  recognized  that  atomization  is  strongly  influenced 
by  upstream  conditions  arising  from  internal  flow 
cavitation,  turbulence,  and  nozzle  geometry.  The 
fidelity  of  the  rate-of- injection  profiles  was  also  shown 
to  be  critical  to  correctly  predict  the  spray  behavior  for 
validation  purposes.  More  recently  and  extending  this 
work,  Senecal  et  al  [6]  has  reported  in  detail  the 
modeling  guidelines  to  simulate  non-evaporating, 
evaporating,  and  reacting  diesel  sprays  utilizing  the 
CONVERGE  solver  in  single-hole  configurations.  In 


this  study,  nominal  ECN  Spray  A  conditions  were  used 
to  demonstrate  key  modeling  features  to  enable  grid 
converged  calculations,  such  as  the  importance  of 
adaptive-mesh-refinement  (AMR),  fixed  grid 
embedding,  and  the  proper  use  of  Lagrangian  parcel 
distribution  for  accurate  hydrodynamic  representation 
of  the  multiphase  flow  [6]. 

Further  studies  demonstrated  the  effect  of  several 
high-resolution  EES  turbulence  techniques  on  resolving 
the  flow  field,  cycle  to  cycle  variations  in  the  context  of 
grid  convergence,  and  its  effect  on  spray  modeling 
constants  [7-8].  Good  agreement  was  found  for  global 
quantities  such  as  liquid  and  vapor  penetration  lengths 
with  experimental  data  for  single-shot  EES  realization. 
However,  deeper  comparison  to  local  quantities  such  as 
velocity  and  mixture  fraction  was  shown  to  require  the 
average  of  many  additional  EES  realizations.  Moreover, 
the  spray  breakup  constants  are  demonstrated  to  also 
depend  on  the  amount  of  resolution  utilized  in  resolving 
the  flow  fields,  hence  the  transition  from  RANS  to  EES 
models  was  reported  to  require  consideration  to 
recalibrate.  Similarly,  Habchi  et  al  [9]  studied  the 
effects  of  EES  ensemble  averaging  using  single-hole 
injectors  in  high-pressure,  high-temperature  vessel  and 
found  good  agreements  with  experimental  data  using 
10-30  EES  realizations  using  the  AVBP  solver.  Hence, 
it  is  generally  agreed  that  using  EES  is  similar  to  a 
single-shot  experimental  injection  requiring  ensemble 
averaging,  while  with  RANS  a  single  shot  realization 
will  be  enough  to  capture  important  spray  quantities 
albeit  with  lower  accuracy. 

The  aim  in  this  study  is  to  extend  the  present 
modeling  scope  to  multi-hole  nozzles  at  the 
standardized  ECN  conditions.  Its  merit  also  includes 
introducing  the  use  of  surrogate  fuel  mixtures  to  model 
evaporating  JP-8  evaporating  sprays  at  diesel 
conditions.  Grid  converged  numerical  calculations  are 
sought  in  order  to  provide  grid  criteria  in  single  and 
multi-hole  nozzle  configurations.  The  models  were 
added  to  the  present  CONVERGE  liquid  fuel  database 
and  validated  extensively  leveraging  the  DOD’s  Spray 
Combustion  Research  Lab  (SCRL)  at  ARL-VTD.  The 
measurements  and  the  experimental  facility  are 
presented  in  more  details  in  the  work  of  Kurman  et  al 
[10]. 

Numerical  Method 

The  CONVERGE  3D-CFD  solver,  developed  by 
Convergent  Science  Inc  (CSI),  has  been  adopted  in  this 
study  to  perform  detailed  spray  simulations  at  realistic 
engine  operating  conditions.  CONVERGE  is  a 
compressible  Navier  Stokes  solver  which  is  based  on  a 
first  order  predictor-corrector  (PISO)  time  integration 
scheme,  and  a  choice  of  second  or  higher  order  finite 
volume  schemes  for  spatial  discretization.  It  features  a 
non-staggered,  collocated,  computation  grid  framework 
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utilizing  a  Rhie-Chow  interpolation  technique  to  avoid 
spurious  oscillations.  An  efficient  geometric  multi-grid 
treatment  is  used  to  solve  the  pressure  equation,  and 
parallel  computing  is  based  on  implementations  of 
either  OpenMP  or  Message  Passing  Interface  (MPI) 
protocols.  It  provides  the  option  of  increasing  resolution 
locally  through  static  fixed-grid  embedding,  and 
dynamically  through  AMR  activated  through  user 
specified  criteria.  Additionally,  it  uses  state-of-the-art 
Eularian-Lagrangian  spray  models,  and  a  parallel 
detailed  chemistry  solver  for  combustion  that  can  be 
fully  coupled  with  CHEMKIN  databases.  The  solver 
also  provides  a  choice  between  a  number  of  modeling 
options  for  the  treatment  of  turbulence,  including  direct 
numerical  simulation  (DNS),  large  eddy  simulation 
(EES)  and  Reynolds-averaged  Navier-Stokes  (RANS). 


Dynamic  Structure  LES  Model  (DSLES) 

In  this  study,  the  gas  phase  is  described  using  the 
Favre  Averaged  Navier  Stokes  Equations  while 
adopting  the  dynamic  structure  LES  model  (DSLES) 
[11-12].  The  compressible  system  of  transport 
equations  for  mass,  momentum  and  energy  are 
presented  here  in  a  LES  framework. 


dp  dpUj 
dt  dXj 


=  0 


dpu, 

dt 


dpupij 

dXj 


dP  da^j  dxij 

dXi  dXj  dXj 


(1) 

(2) 


where  the  subgrid  stress  tensor  is  given  by 

Tjy  =  p(u^j  —  ufij)  (3) 


In  the  DSLES  methodology,  the  subgrid  stress 
tensor  is  modeled  through  a  non-viscosity  based  one- 
equation  LES  model.  The  transport  equation  for  the 
subgrid  kinetic  energy  is  written  as. 


dpk  dpUjk 
dt  dXj 


dXj 


p  dk\  _ 


Pr  dx- 

^  *  sgs 


T  -  S 


■pe 
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where  the  subgrid  kinetic  energy  is  given  by. 


k  =  —  {upij  —  ufij) 


(5) 


and  the  subgrid  dissipation  rate  is  given  by. 


(6) 


Model  constants  in  this  formulation  is  Q  =  1,  A 
represents  the  local  grid  length  scale,  S^J  is  the  filtered 
strain  rate  and  Pr^gs  is  set  1.39.  The  subgrid  stress 
tensor  is  now  modeled  as. 


Ti;  =  2kp-^  (7) 

where  the  Leonards  stress  tensor  is  classically 
defined  as,  Ljy  =  (ujiZ]  —  rqu])  including  the  test  filter 
and  grid  filter  level  in  its  definition. 

DSLES  methodology  requires  the  scalar  and 
energy  transport  equations  be  written  using  the  eddy 
viscosity  concept. 


dt  dXj 


d  (  dY, 

pW  +  d,) 


dXj 


dXj 


(8) 


where  —  C^k°'^  /Sc  defines  the  turbulent 
diffusion  parameter  with  Q  =  0.5  and  Schmidt  number 


5c  =  0.78. 


Time  Advancement 

Equations  (1)  and  (2)  are  solved  numerically  using 
a  classical  predictor-corrector  scheme  in  which  the 
velocity  field  is  first  integrated  using  the  Navier  Stokes 
equation  (2)  and  then  corrected  to  enforce  mass 
conservation  (1)  using  a  modified  pressure.  This  is 
achieved  using  the  modified  Pressure  Implicit  with 
Splitting  of  Operator  (PISO)  algorithm  first  introduced 
by  Issa  [13]. 

An  illustration  of  the  PISO  predictor-corrector 
sequence  is  presented  below  (note  that  H*  represents 
convection,  diffusion,  S  accounts  for  imp licit/exp licit 
sources,  and  denotes  the  intermediate  steps)  for  a  one¬ 
time  correction  procedure. 


Predictor  step: 

p"u*  p"uf  _  dP" 
dt  dt  dxi 

Corrector  step: 

p*u**  p^uf  _  dP* 
dt  dt  dxi 


(9) 


(10) 


d^ 

dxidxi 


(P*  -  P")  - 


(P*  -  P^)(p 
dt^ 


dp'^u*  _ 
dxi  )  dt 


Update  step: 


p*ul*  p^u-l 
dt  dt 
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dxi 


(P*  -  P") 


(11) 
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where  dt  is  the  numerical  integration  time  step. 
Once  the  pressure  is  solved,  the  velocity  is  re-calculated 
according  to  the  updated  momentum  equation. 

Physical  Spray  Models 

Spray  modeling  is  treated  using  the  “blob” 
injection  method  of  Reitz  and  Diwakar  [14].  Blobs  of  a 
characteristic  size  are  injected  following  a  statistical 
distribution  into  the  computational  domain.  Primary 
and  secondary  breakups  are  subsequently  simulated 
based  on  the  Kelvin-Helmholtz  (KH)  and  Rayleigh- 
Taylor  (RT)  instability  methods.  Note  that  the  breakup 
length  is  not  determined  a-priori  (breakup  length 
concept)  and  is  calculated  as  a  part  of  the  solution. 

In  the  KH  wave  model,  atomization  is  treated  using 
stability  analysis  for  liquid  jets.  The  breakup  of  the 
injected  blobs  and  resulting  drops  of  radius  is 
calculated  by  assuming  that  the  drop  radius  is 
proportional  to  the  wavelength  of  the  fastest  growing 
unstable  surface  wave  A^fj. 

It  is  written  as,  r  =  where  Bg  is  a  model 

constant.  Droplet  size,  and  its  rate  of  change  is 
calculated  as, 

drp  _  (tq  -  r) 

dt  Tjffj  (12) 

where  the  breakup  time  constant,  ,  is  calculated  as. 


3.726Rj^ro 
^KH  Hkh 


(13) 


and  the  maximum  growth  rates  and 

corresponding  wavelengths  A^fj  have  been  simplified 
and  defined  as  follows, 

fpia^\  0.34 -I- O.SSfRe^-^ 
j  ^  (1  -f  Z)(l  -f 

Akh  (1  +  0ASZ°-^)(1  +  0AT°'^) 

=  9  02  - _ — - _ - 

“  '  (l -f  0.871Rei-67)°  ®  (15) 

where, 

Z  =  We^-^/Rei ,  T  =  We°-^,  VRe,  =  piU^aja , 
WCg  =  pgU^a/a,  and  Re;  =  Bajvi 

The  present  Rayleigh-Taylor  mechanism 
formulation  includes  viscosity  variations  in  the  growth 
rate  equation. 


where  knj  is  the  wavenumber,  pi  is  the  liquid 
viscosity,  Pg  is  the  gas  viscosity,  p;  is  the  liquid 
density,  Pg  the  gas  density,  a  is  the  deceleration  of  the 
drop,  and  a  is  the  liquid  surface  tension.  The  wave 
number  corresponding  to  the  maximum  growth  rate 
Kjfi-  —  2nlA[fY  is  calculated  through  a  bisection  method 
with  equation  (16).  The  value  is  updated  to  calculate  the 
maximum  growth  rate  D-ht  .  The  predicted  RT  model 
drop  size  is  then  expressed  as, 

fRT  —  ^RT^KH  (12) 

where  is  the  model  constant,  and  A^g  is  the 
predicted  RT  wavelength.  In  summary  the  spray  model 
constants  used  in  this  study  are  defined  as:  Bp  =  0.6 , 
Ri  =  10  ,  Cuf  =  0.1 ,  and  Q  =  1. 

Lastly,  other  spray  processes  that  were  modeled 
include  droplet  distortion  and  drag,  droplet  interactions 
in  terms  of  collision  and  coalescence,  turbulent 
dispersion,  and  evaporation.  The  two  phases  are 
coupled  through  the  exchange  of  mass,  momentum,  and 
energy,  represented  in  the  appropriate  source  terms  in 
the  gas-phase  equations.  These  models  are  reviewed 
and  described  in  detail  in  a  recent  report  by  Bravo  et  al 
and  the  reader  is  referred  to  that  article  [15]. 

Experimental  Data  and  Compntational  Setnp 

Model  validation  is  carried  out  with  our  in-house 
database  at  ARL-VTD,  Spray  Combustion  Research 
Laboratory  (SCRL).  This  is  a  state-of-the-art  facility 
housing  a  HTPV  providing  the  high  temperature  and 
high  pressure  thermodynamic  conditions  obtained  in  a 
diesel  engine.  Unlike  CVP  (Constant  Volume  Prebum) 
vessels  this  rig  does  not  require  a  pre-burn  phase  to 
obtain  this  condition  and  rather  it  is  comprised  of  four 
subsystems  including:  gas  compressor,  gas  heater,  test 
vessel,  and  control  system  to  achieve  nearly  quiescent 
and  steady  thermodynamic  conditions.  An  onsite 
nitrogen  generator  (Model  N-350  PSA)  and  a  high- 
pressure  compressor  (Sauer  Model  WP4351) 
continuously  deliver  high-pressure  nitrogen  to  the 
HTPV.  A  2-stage  cartridge  electrical  heater  and  a 
ceramic  heater  are  used  to  increase  the  gas  temperature 
with  the  use  of  a  control  system  (closed  loop  PID)  to 
adjust  for  both  the  chamber  pressure  and  the  power  of 
the  heater  to  obtain  homogenous  temperature 
conditions.  To  characterize  the  liquid  and  vapor  global 
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measurements  an  advaneed  optieal  set  up  with  an  LED 
array  is  used  for  Mie  scattering  technique  and  single 
LED  lighting  is  used  for  schlieren  technique.  Bosch 
common-rail  fuel  injectors  with  147  jxm  nominal 
diameter  were  used  in  experiments.  X-ray  CT  scan 
images  of  the  injectors  are  shown  in  Figurel. 


Parameter 

Quantity 

Injector  Types 

BOSCH  CRIN3 

Number  of  Holes 

1 

2 

Nominal  Orifice 
Diameter  (/rm) 

147 

147 

Fill  Gas 

N2  99% 

Chamber  Density  (kg/m^*) 

22.8 

Chamber  pressure  (bar) 

60 

Fuel 

JP-8 

Fuel  Temperature  (K) 

373 

Fuel  injection  Quantity 
(m.g) 

7.05 

11.52 

Injection  Pressure  (bar) 

1000 

Injection  Duration  (ms) 

1.1 

Pulse  Duration  (ms) 

0.7 

Table  1.  Test  matrix  combustion  vessel  and  fuel 
injector  with  experimental  conditions  for  CFD  model 
validation. 


Figure  1.  Rendered  X-ray  CT  scans  (top  row)  with 
corresponding  image  slices  (bottom  row)  for  single 
(a)  and  two-hole  (b)  injectors. 

It  is  noted  that  the  single-hole  nozzle  orifice 
presented  in  Figure  la  has  been  manufactured  using 
wire-EDM  technique  while  welding  shut  the  side 
orifices.  The  top  row  is  a  rendered  image  from  the  CT 
scan  showing  that  the  orifice  is  not  aligned  with  axial 
position.  Two-hole  nozzle  orifices  also  have  welding  on 
the  side  orifices  to  allow  control  of  the  fuel  injection 
orifices.  The  nominal  distance  between  orifices  in  this 
case  is  approximately  1  mm  with  60  degree  internal 
angle  as  shown  in  the  Figure  lb.  The  images  in  Figure  1 


are  insightful  in  presenting  manufacturing  eccentricities 
that  are  useful  for  full  undersfanding  of  the  spray 
behavior. 

To  characterize  the  rate  of  injection  (ROI)  and  the 
injected  fuel  mass,  each  fuel  injector  was  tesfed  in  an 
lAV  fuel  injection  analyzer.  For  the  injector 
characterization  presented  in  this  work,  a  100  shot 
average  was  acquired  for  the  injected  fuel  mass,  ROI, 
and  current  profile.  The  ROI  profiles  used  in  the  CFD 
work  and  their  comparison  to  experiments  are  presented 
in  Figure  2.  Other  operating  conditions  used  in  this 
study  are  provided  in  Table  1. 


Figure  2.  Experimental  Rate  of  Injection  Profiles 
compared  with  the  profiles  used  in  CFD  work. 


In  this  study  we  use  numerical  simulations  to 
model  transient  non-reacting  JP-8  sprays  under  diesel 
engine  conditions.  The  CONVERGE  multiphase  solver 
is  used  invoking  a  classical  Lagrangian/Eularian 
formulation  to  treat  the  gas-phase  and  liquid  phase  with 
wave-based  spray  models.  The  study  includes  single¬ 
hole  and  two-hole  nozzle  configurations  at  the  nominal 
Spray  A  ambient  conditions. 

There  are  several  important  procedures  to  be 
mindful  when  running  spray  simulations.  First,  an 
AMR  procedure  must  be  adopted  to  control  the 
refinement  levels  in  the  vicinity  of  the  spray.  Figure  3 
shows  the  effects  of  AMR  refinements  in  the  near  the 
spray  calculated  with  a  velocity  criteria  using  a  base 
grid  of  2  mm.  In  this  study,  AMR  is  activated  for  the 
velocity  field  and  its  use  is  critical  in  keeping  the  cell 
counts  at  realistic  operating  levels.  Note  also  that  fixed 
grid  embedding  (in  spray  liquid  core  region)  is  used  to 
define  the  reference  minimum  cell  sizes,  dx,  and  the 
embedding  scale.  The  parameters  are  related  through 
the  base  grid  cell  size,  dXf,ase^  follows:  dx  = 
The  coarsest  base  grid  cell  size 
of  2  mm,  and  an  embed  scale  of  3  corresponds  to 
dx^O.25  mm.  When  conducting  grid  refinements,  the 
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number  of  injected  particles  must  be  increased  at  each 
refinement  level.  Using  a  fixed  number  of  injected 
particles  will  lead  to  an  artificial  increase  of  the  liquid 
penetration  length  as  was  reported  in  a  previous  study 
by  Senecal  et  al  [6].  This  is  due  to  the  decrease  of 
relative  gas  to  liquid  mass  that  occurs  in  each  cell  as  the 
mesh  is  refined.  Overlooking  this  can  artificially  cause 
the  gas  velocity  to  approach  the  particle  velocity  at 
every  cell  leading  to  an  artificial  decrease  in  particle 
drag. 


2  mm 


108  mm 


Figure  3.  Computational  mesh  showing  the  effects  of 
AMR  near  the  spray  region,  and  a  coarse  base  mesh 
outside  of  it. 

Table  2  documents  the  number  of  resolution  levels, 
denoted  as  embedding  scale,  and  the  critical  number  of 
injected  particles  used  in  grid-sensitivity  analysis. 
Embed  scales  3-5  were  used  in  this  study. 


Embed 

Scale 

Cell  Size 
dx(mm) 

d„/dXi 

Injected 

Parcels 

2 

0.5 

0.294 

12,500 

3 

0.25 

0.588 

50,000 

4 

0.125 

1.176 

200,000 

5 

0.0625 

2.352 

800,000 

Table  2.  AMR  embed  scale  (column  1),  cell  sizes 
(column  2),  number  of  cells  inside  of  the  nozzle 
diameter  (column  3),  and  number  of  injected  parcels 
(column  4)  used  for  evaporating  sprays. 


The  guidelines  above  were  utilized  in  conducting 
grid  convergence  studies  to  demonstrate  the  numerical 
accuracy  of  the  model.  Numerical  convergence  can  be 
achieved  at  cell  sizes  0.125  mm  (1  hole)  and  0.0625 
mm  (2  hole)  configurations.  Physical  validation  is 
conducted  by  comparing  grid-converged  results  to  the 
measurements.  In  this  study  both  liquid  and  vapor 


penetration  data  was  available  for  this  task  as  well  as 
animations  showing  plume  structure. 

Overview  and  Selection  of  JP-8  Surrogate  Fuel 

There  exist  several  categories  of  aviation  fuels  in 
operation  ranging  from  commercial  aviation  fuels, 
including  Jet  A  and  Jet  A-1,  to  military  jet  fuels  using  a 
different  classification  system  JP  numbers  spanning 
from  JP-1  to  JP-8  varying  with  application  [16].  A 
general  summary  of  fuel  property  data  is  given  by 
Edwards  et  al.  in  his  review  of  relevant  military  fuels 
[17].  Fuel  JP  4  is  reported  as  a  gasoline/kerosene  blend 
used  by  the  U.S.  Air  Force  (USAF)  until  the  1980’s 
when  it  was  replaced  by  JP-8;  a  kerosene  based  fuel. 
JP-8’s  specifications  are  similar  to  Jet  A  (Jet  A-1), 
which  is  a  commercial  aviation  kerosene  fuel,  with  the 
difference  due  to  the  addition  of  three  additives;  the 
additives  are:  a  lubricity  improver/corrosion  inhibitor, 
antistatic  additive,  and  a  fuel  system  icing  inhibitor. 
Fuel  JP-5  is  a  high  flash-point  kerosene  used  for  aircraft 
flying  from  U.S.  Navy  ships,  while  fuel  JP-7  is  an 
exclusive  kerosene  fuel  developed  for  supersonic 
applications  related  to  the  military  SR-71  “Blackbird” 
aircraft.  The  following  table  lists  average  property  data 
for  the  aforementioned  military  fuels  [17]. 


Property 

JP-4 

JP-7 

JP-8 

Approx  formula 

^8.5  ^17 

^12^25 

^11^21 

H/C  ratio 

1.99 

2.02 

1.91 

Boil  Range,  F 

140-460 

370- 

480 

330-510 

Freeze  Point,  F 

-80 

-47 

-60  JP8/JetAl 
-50  Jet  A 

Flash  point,  F 

-10 

140 

127 

Net  heating  val 
(BTU/lb) 

18,700 

18,530 

18,580 

Specific  gravity 
60  F 

0.76 

0.79 

0.81 

Critical  T,  F 

450 

305 

340 

Critical  P,  psia 

450 

305 

340 

Average  comp. 

aromatics,  vol% 

10 

19 

18 

naphthenes 

29 

32 

20 

paraffins 

59 

65 

60 

olefins 

2 

- 

2 

sulfur,  ppm 

370 

60 

490 

Table  3.  Bulk  fuel  properties  of  typical  aviation 
military  fuels  (Edwards  et  al  [16]). 


An  inspection  of  Table  2  reveals  that  the  major 
hydro-carbon  constituents  found  in  JP-8  are:  paraffins, 
60  vol%,  followed  by  naphthenes  (cycloparaffins),  20 
vol%.  It  is  reported  that  paraffins  are  composed  of 
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primarily  of  straight-chained  or  branched  molecules, 
with  mono-cycloparaffins  also  being  significant.  A 
major  paraffin  component  is  usually  n-decane  or  n- 
dodecane,  characterized  with  high  H/C  ratios  and  high 
energy  release  rates  per  unit  weight.  Naphthenes  have 
lower  H/C  ratios  and  exhibit  lower  gravimetric  heating 
values;  however,  they  feature  higher  volumetric  energy 
densities  and  lower  freezing  points  providing 
operational  advantages.  Aromatic  content  in  JP-8  fuel  is 
of  18  vol%  and  constitutes  the  third  largest  class  of 
compounds.  Aromatic  hydrocarbons  can  be  monocyclic 
(MAH)  and  polycyclic  (PAH),  where  polycyclic- 
aromatic  hydrocarbons  are  well-known  soot  precursors 
and  is  an  undesirable  result  of  incomplete  combustion 
directly  affecting  engine  performance.  Government 
regulations  for  aviation  fuels  limit  the  aromatic  content 
in  JP-8  to  be  around  18%  -  25%.  Lastly,  trace  species 
such  as  sulfur  and  metals  is  important  in  fuels  for 
thermal  stability  and  emissions.  The  average  sulfur 
content  among  JP-8  fuels  is  tabulated  as  490  ppm 
(0.049  mass%),  less  than  the  0.3  mass%  government 
mandated  maximum. 


Surrogate-C  JVfj  Xj 

n-dodecane 

C12H26 

170 

0.6 

0.714 

Methylcy¬ 

clohexane 

C7H14 

98 

0.2 

0.137 

o-xylene 

0 

0 

106 

0.2 

0.148 

H/C  ratio 

1.97 

^mix 

142.8 

Surrogate-Mod 

Aachen 

Mi 

Xi 

Vi 

n-dodecane 

C12H26 

170 

0.8 

0.85 

trimethyl- 

benzene 

C9H12 

120 

0.2 

0.15 

H/C  ratio 

1.99 

^mix 

160 

_  1 

Table  4.  Table  4.  JP-8  surrogate  composition.  Chemical 
formula  (column  2),  molecular  weight  (column  3), 
mole  fraction  (column  4),  mass  fraction  (column  5). 
Note  the  last  rows  provide  each  surrogate  H/C  ratio  and 
molecular  weight. 

In  this  study,  the  use  of  surrogate  fuels  is  adopted 
in  an  effort  to  model  JP-8  liquid  physical  properties.  A 
literature  review  reveals  that  there  is  a  need  for  a 
complete  characterization  of  JP-8  liquid  phase 


properties  (i.e.,  vapor  pressure)  that  spans  an  extensive 
temperature  range  for  use  in  CFD  models.  The  models 
used  in  this  study  were  selected  based  on  literature 
review  on  gas-phase  JP-8  surrogates  and  their 
availability  [16].  They  are  a  two-component  Modified 
Aachen  surrogate  with  constituents  of  n-dodecane  and 
trimethylbenze,  and  a  three-component  Surrogate-C 
with  constituents  of  n-dodecane,  o-xylene,  and 
methylcyclohexane.  The  exact  mole  and  mass-fraction 
compositions,  mixture  molecular  weight  and  hydrogen 
to  carbon  ratios  (H/C)  are  presented  in  Table  4. 

Results  and  Discussion 

In  this  section  we  present  modeling  results  of  JP-8 
evaporating  spray  simulations  with  single-hole  and  two- 
hole  nozzle  configurations  at  the  Spray  A  ambient 
conditions.  Using  a  RANS  approach  we  compare  the 
spray  global  parameters,  liquid  and  vapor  penetration 
profiles,  and  study  the  effect  of  resolution  and  surrogate 
fuels  on  spray  structure.  With  the  LES  approach  we 
investigate  further  the  turbulent  flow  field  mixing 
behavior  surrounding  the  atomizing  jet  through 
visualization  of  velocity  profiles,  vorticity,  and  mixture 
fraction. 

For  reference,  the  calculations  of  liquid  and  vapor 
penetrations  for  simulations  are  defined  as  follows: 


Parameter 

Definition 

Liquid 

Penetration 

Axial  distance  encompassing  97%  of  the 
injected  liquid  fuel  mass. 

Vapor 

Penetration 

Maximum  axial  distance  from  the 
injector  orifice  where  the  fuel  mass 
fraction  Yf  =  0.1%. 

Table  5.  Definitions  of  current  evaporating  spray  global 
parameters. 


Surrogate  Grid  Convergence  (Single-Hole  Nozzle) 

Of  interest  here  is  to  develop  a  model  in  which  the 
global  diagnostics  presented  above  asymptotically 
converge  as  the  mesh  is  refined.  Grid  convergence  in 
this  study  has  been  achieved  by  utilizing  the  solver’s 
ability  to  refine  the  grid  dynamically  by  using  AMR, 
and  locally  by  use  of  fixed  grid  embedding  modules.  A 
correct  representation  of  the  particle  hydrodynamics  is 
also  important  in  predicting  the  spray  behavior.  This  is 
achieved  by  increasing  the  number  of  injected  particles 
with  resolution  level.  The  convergence  grid  parameters 
selected  here  correspond  to  grid  cell  sizes  of  dx  — 
0.25, 0.125, 0.0625  mm.  In  this  case,  the  relative 
nozzle-to-cell  ratio  increases  approximately  by  two 
with  each  level  reaching  a  peak  value  of  = 

2.352  for  dx  —  0.0625  mm  . 

The  first  set  of  grid  sensitivity  results  are 
conducted  with  Surrogate-C  showing  good  comparison 
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to  measurements  for  liquid  and  vapor  penetration 
profiles  as  seen  in  Figure  4a-b.  In  Figure  4a  (top),  for  a 
cell  size  dx  —  0.25  mm,  there  is  a  sharp  peak  in  the 
transient  region  of  the  liquid  penetration  profiles.  This 
unphysical  spike  in  penetration  is  a  clear  indicator  of 
unresolved  grid  effects.  Improvements  are  gradually 
obtained  for  cells  =  0.125,  0.0625  mm  where  the 
finer  grid  removes  the  peak  and  hence  is  able  to  capture 
the  transient  spray  behavior  more  accurately.  Surrogate- 
C  vapor  penetration  behavior  also  shows  asymptotically 
convergence  behavior  as  resolution  is  increased,  see 
Figure  4b  (bottom).  In  this  case,  grid  effects  become 
evident  after  a  time  of  0.2ms.  Some  of  the  reported 
discrepancies  with  the  present  vapor  phase  length  are 
attributed  to  the  uncertainties  in  the  numerical  length 
calculation. 


Figure  4.  Surrogate-C,  single  hole  nozzle,  grid 
convergence  behavior  with  liquid  (a)  and  vapor  (b) 
penetration  profiles. 

A  further  study  utilizing  the  Modified  Aachen 
surrogate  shows  similar  results  as  with  Surrogate-C  in 
terms  of  the  liquid  and  vapor  profiles.  Figures  5a-b.  The 
similarities  in  physical  behavior  for  each  surrogate  can 


be  attributed  to  its  small  variance  (less  than  2%)  while 
both  values  are  near  the  IP-8  reported  value ///C  = 
1.99.  It  is  also  important  to  note  that  the  temporal 
breakup  constant,  Bi ,  was  recalibrated  to  represent  the 
present  spray  behavior  with  IP- 8  (repeated  here  for 
clarity)  =  3.726Biro/A^fjD.^YiThis  was  achieved 
by  carrying  out  a  parametric  study  with  and 
selecting  the  best  approximation  with  experiments.  In 
this  process  the  breakup  constant  was  modified  for  JP-8 
surrogates  from  B-^  —  7  (diesel  fuels)  to  =  10  . 


Figure  5.  Surrogate-Modified  Aachen,  single  hole 
nozzle,  grid  convergence  behavior  with  liquid  (a)  and 
vapor  (b)  penetration  profiles. 

Figure  6  presents  the  computed  spray  structure  at  times 
50,  100,  and  200  jis  after  start  of  injection  at  a 
resolution  level  dx  =  0.0625.  The  vertical  line  marks 
the  computed  (PANS)  liquid  lengths.  Note  that  the 
PANS  turbulence  model  adopted  in  this  study  captures 
very  well  the  transient  and  steady  penetration  behavior; 
however,  higher  fidelity  is  needed  to  capture  the  mixing 
spray  behavior  with  more  accuracy. 
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Figure  6.  Single  hole  nozzle,  Lagrangian  particles 
depicting  penetration  layer  growth  using  the 
Modified  Aachen  Surrogate  at  c&  =  0.0625  mm 


Surrogate  Grid  Convergence  (Two-Hole  Nozzle) 

The  effect  of  multi-hole,  in  this  case  two-hole 
adjacent  nozzle,  is  investigated  to  numerically  seek 
grid-convergence  behavior.  A  grid  converged  model 
will  define  the  resolution  required  in  carrying  out  well 
resolved  simulations.  This  configuration  features  two 
emerging  plumes  at  different  orientations.  They  are 
denoted  in  terms  of  degrees  as  plume  0  degree  (vertical) 
and  plume  60  degree  (adjacent),and  both  have  a 
manufacturer  provided  included  angle  of  159  degrees. 

Figure  7a-b  presents  the  effect  of  spray  liquid  and 
vapor  behavior  with  increased  resolution  and 
comparison  to  experiments  using  Surrogate-C.  The 
resolution  levels  selected  are  for  dx  =  0.25,  0.125,  and 
0.0625  mm.  Figure  7a  (top)  shows  a  variations  in  liquid 
phase  length  for  the  0  deg  and  60  deg  plumes  with 
coarse  resolution.  This  is  an  effect  of  grid  to  spray 
misalignment  which  improves  as  the  grid  is  refined.  For 
a  cell  size  of  0.0625  mm  both  plumes  show  similar 
behavior  in  liquid  phase  penetration  establishing  the 
grid  criteria  for  two-hole  plumes.  At  the  finest  scales. 
Figures  7a-b,  both  show  differences  with  the 
measurements.  The  discrepancies  in  measurements  and 
modeling  are  for  2  mm  and  4  mm  in  the  liquid  phase 
and  vapor  phase,  respectively.  Note  that  this  is  the 
current  results  using  the  RANS  methodology  and  will 
be  interrogated  further  with  LES;  however,  RANS  is 
able  to  capture  the  overall  behavior  with  a  good  degree 
of  accuracy. 

Figure  8  shows  the  behavior  of  the  Modified 
Aachen  surrogate  from  the  same  grid  sensitivity  study. 
Similar  results  are  obtained  showing  large  discrepancies 
for  coarse  cell  (dx  =  0.25  mm)  sizes,  improvements  for 
dx  =  0.125  on  the  60  degree  plume,  and  a  complete 
attenuation  of  the  grid  effect  for  the  finest  cell  size 


Figure  7.  Surrogate-C,  two-hole  nozzle,  grid 
convergence  behavior  with  liquid  (a)  and  vapor  (b) 
penetration  profiles. 


#  Measurements  (Plume  Odeg) 

—A -  Measurements  (Plume  60deg) 

.  dx  =  0.25  mm,  Plume  0  deg,  Surrogate  Mod  Aachen 

-  dxs0.25mm.  Plume  60  deg,  Surrogate  Mod  Aachen 

dx  =  0.125  mm,  Plume  0  deg,  Surrogate  Mod  Aachen 

-  dxs  0.125  mm,  Plume  60  deg.  Surrogate  Mod  Aachen 

.  dx  s  0.0625  mm.  Plume  0  deg.  Surrogate  Mod  Aachen 

-  dx  =  0.0625  mm,  Plume  60  deg.  Surrogate  Mod  Aachen 


ASOI  (ms) 


Figure  8.  Surrogate-Modifed  Aachen,  two  hole 
nozzle,  grid  convergence  behavior  with  liquid 
penetration  profile 


(dx  =  0.0625  mm).  Vapor  phase  grid-analysis  results 
show  similar  behavior  as  with  Surrogate-C. 

Figure  9  presents  the  computed  spray  structure  at 
times  5,  12,  and  40  fis  after  start  of  injection.  The 
liquid  phase  behavior  is  clearly  demonstrated  for  each 
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nozzle  orifice  with  the  correct  orientation.  The  vertical 
line  drawn  on  each  plume  marks  the  computed  liquid 
lengths  for  each  case.  Note  that  this  is  the  result  of  a 
RANS  turbulence  modeling  approach  with  Lagrangian 
particle  spray  visualization. 

9  mm 


15  mm 


Figure  9.  Surrogate-Modified  Aachen,  two  hole 
adjacent  orifice  nozzle,  Lagrangian  particles 
depicting  liquid  penetration  layer  growth  in  spray 
formation,  dx  =  0.0625  mm. 

Large  Eddy  Simulations  Results 

This  section  present  the  results  obtained  using  the 
dynamic  structure  LES  model  described  in  the 
numerical  methodology  section.  A  fine  base  grid  cell 
size  of  dx  -  0.0625  mm  with  800,000  total  injected 
parcels  for  each  plume  was  utilized,  AMR  criteria  for 
evaporating  sprays  was  based  on  velocity  variations.  A 
total  of  9  independent  simulations  (9  cycles)  were 
carried  out  for  each  experiment  resulting  in  averaged 
profiles  for  liquid  and  vapor  phase  penetrations.  This 
choice  was  based  on  the  previous  work  by  Senecal  et  al 
[7]  where  it  was  reported  that  a  total  of  9  cycles  were 
needed  as  the  minimum  number  of  LES  realizations  to 


maintain  a  reasonable  representation  of  the  average 
spray  behavior  for  the  Spray  A  configuration  (the 
average  of  9  LES  realization  will  have  a  mean  absolute 
error  of  less  than  +!-  2.0  m/s  and  0.005  for  velocity  and 
mixture  fraction,  respectively,  95%  of  the  time).  In  this 
section,  JP-8  fuel  is  modeled  by  selection  of  the  Aachen 
surrogate.  Note  also  the  LES  spray  breakup  time 
constant  used  was  Bi  —  5. 

Single-hole  LES  results  are  presented  in  Figure 
lOa-b  showing  each  LES  cycle  realization  in  light  red 
line  along  with  the  averaged  quantities  in  solid  red  line 
and  comparison  with  experiment. 


Figure  10.  Single  hole  LES  results  using  the  Dynamic 
Sti-ucture  model  (dx  =  0.0625  mm)  liquid  phase  (a)  vapor 
phase  (b)  penetration 

For  the  single-hole  nozzle  the  results  show  good 
agreements  for  liquid  phase  length  profiles,  and 
moderate  success  with  vapor  phase.  The  vapor  phase 
discrepancies  are  larger  at  earlier  times  and  decrease  to 
less  than  2%  difference  at  1  ms  (this  may  be  due  to  the 
use  of  single  component  fuel  compared  to  the  real  fuel 
in  the  experiment  (Different  evaporation  rates).  At  this 
point  the  modeling  results  are  within  experimental 
uncertainty.  Additional  single-hole  results  are  presented 
and  visually  compared  to  measurements  in  Figure  11a- 
b.  The  experimental  data  shows  the  liquid  and  vapor 
phase  lengths  by  using  a  Mie  (liquid)  and  schlieren 
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(vapor)  technique  for  visualization  and  are  shown  side 
by  side  with  LES  results.  In  this  Figure,  the  modeling 
results  present  the  velocity  contours  and  shows  a  peak 
velocity  of  500  m/s  at  the  injector  nozzle  followed  by  a 
typical  radial  Gaussian  velocity  distribution  and  decay 
in  the  axial  direction.  The  vapor  phase  here  is  depicted 
as  a  solid  black  line  surrounding  the  turbulent  plume 
region,  the  imaging  times  are  33,  165,  495,  891,  1188 
[IS.  Figure  12a-d  presents  the  spray’s  overall  mixing 
and  scalar  behavior  in  terms  of  the  mixture  fraction, 
temperature,  turbulent  kinetic  energy,  and  vorticity  at 
the  same  time  intervals  as  Figure  11.  In  Figure  12a, 
mixture  fraction  distributions  also  follow  a  Gaussian 
type  profde  with  peak  values  (Z  =  0.4)  at  the  centerline 
near  the  orifice  and  decaying  radially  at  fixed  times  and 
axially  as  time  progresses.  Figure  12b  shows  the  pre¬ 
heated  fuel  entering  the  chamber  with  a  rail  pressure  of 
1000  bar  while  undergoing  strong  mixing  with  heated 
ambient  at  Too  =  900  K  .  To  elucidate  further  on  the 
turbulence  transport  and  mixing  process.  Figure  12c-d 
show  the  turbulent  kinetic  energy  and  vorticity 
contours.  The  contours  show  the  turbulent  kinetic 
energy  mechanism  to  be  dominant  in  the  spray 
centerline  due  to  the  high  pressure  injection  process. 
The  emerging  spray  interacts  with  quiescent  ambient  air 
which  in  turn  generates  the  3D  surface  instabilities  that 
lead  to  the  spray  breakup  behaviors.  Vorticity 
dominates  in  this  core  region  due  to  the  entrainment  of 
air,  and  mixing,  as  a  result  of  the  injected  fuel.  This 
analysis  provides  a  complete  picture  of  the  spray 
breakup  and  mixing  process  for  single-hole  nozzle  case. 

The  two-hole  nozzle  orifice  was  also  modeled 
using  the  LES  approach.  Equal  spatial  resolution  (dx  = 
0.0625  mm)  and  time  averaging  (9  cycles)  techniques 
were  utilized  as  the  single-hole  case  to  provide 
comparable  results.  The  number  of  injected  parcels  was 
kept  fixed  as  the  single-hole  value  of  800,000  per  hole 
or  total  1.6  million  in  order  to  accurately  resolve  the 
parcel’s  hydrodynamic  effects  on  the  mesh.  In 
particular,  we  note  from  experiments  that  scaling  of  the 
injected  mass  is  non-linear  leading  to  a  relative 
decrease  in  rate  of  injection  (ROI)  peak  profile 
magnitude.  It  is  also  observed  that  the  ROI  profiles 
used  in  this  study  was  the  average  from  the  two  orifices 
as  the  present  measurement  system  provides  total 
quantities. 

Figure  13a-b  shows  the  instantaneous  liquid 
penetration  lengths  in  red  light  lines,  alongside  the 
mean  profile  in  red  solid  line,  and  its  comparison  to 
uncertainty  quantified  measurements.  Note  that  the 
modeled  instantaneous  penetration  profiles  feature 
larger  fluctuations  in  liquid  phase  length  when 
compared  to  the  single-hole  profiles  for  both  plumes. 
This  could  arise  from  potential  insufficient  parcel 
number  distributions  not  scaling  with  the  values  from 


the  single-hole  case.  The  averaged  value  is  consistent 
between  both  plumes  at  20  mm  and  is  an  over 
prediction  to  experiments  by  3-4  mm  to  the  uncertainty 
region.  When  compared  to  single-hole  experiments,  the 
mean  LES  liquid  penetration  profiles  show  similar 
behavior.  This  implies  that  there  is  not  much  of  an 
impact  in  liquid  length  when  injecting  with  single  and 
multi-hole  nozzles  at  this  particular  orientation  (60 
degrees,  approximately  1  mm  apart). 


Figure  13.  Two  hole  nozzle,  LES  liquid  penetration 
profiles  using  the  Dynamic  Structure  model  (dx  =  0.0625 
mm)  for  plume  0  deg  (a)  and  plume  60  deg  (b). 

Figure  14  shows  an  ensemble  of  9  vapor  length 
LES  realizations,  the  resulting  LES  mean  profile,  and 
its  comparison  to  experiments.  There  are  significant 
differences  (~  8  mm)  in  the  vapor  mixing  behavior  with 
experiments  which  results  in  a  similar  trend  observed  in 
two-hole  RANS  simulations.  In  Figure  15  the 
calculations  are  compared  to  experimental  Mie 
Scattering  images  in  an  axial  orientation  at  selected 
times. 
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Figure  14.  Two  hole  nozzle,  LES  vapor  penetration 
profiles  using  the  Dynamic  Structure  model  {dx  = 
0.0625  mm) 

The  interpretation  of  the  results  between  single¬ 
hole  and  multi-hole  nozzle  configuration  requires 
consideration  of  several  critical  injector  components. 
Atomization  in  sprays  is  highly  dependent  on  the 
upstream  conditions  that  occur  inside  the  nozzle 
structure  such  as  cavitation  and  turbulence  (physical 
effects)  and  fluid-structure  interaction  (interfacial 
effect).  The  internal  nozzle  flow  and  passage  through 
orifices  plays  a  critical  role  in  determining  the  rate  of 
injection  profiles,  injected  mass,  and  other  relevant 
spray  quantities.  Needle  motion  in  the  lateral  direction 
(wobble)  has  been  reported  for  rail  pressures  of  1000 
bar  or  greater.  This  effect  has  a  strong  impact  on  the 
emerging  spray  as  the  needle  wobble  can  potentially 
and  cyclically  obstruct  fuel  passageways.  For  sprays 
emerging  from  two  adjacent  orifices,  there  is  an 
additional  geometrical  asymmetry  leading  to  pressure 
disturbances  that  may  be  preventing  the  flow  to  be 
distributed  into  the  orifice  channels  uniformly.  In 
Lagrangian  modeling  of  sprays,  as  is  the  present  effort, 
the  spray  is  specified  from  point  sources  and  imposing 
a  Bernoulli  type  injection  velocity  derived  from  ROI 
profiles,  injection  quantities,  and  discharge 
characteristics.  Hence,  the  internal  flow  details 
occurring  upstream  of  the  spray  injection  region  are  not 
modeled.  It  is  thought  that  these  effects  are  prominent 
for  two-hole  (asymmetrical)  nozzles  where  the 
interaction  between  the  needle  and  displaced  fluid  and 
orifice  location  is  vastly  different  than  the  single-hole 
case. 

Conclusions 

In  this  study  we  have  conducted  validated 
numerical  simulations  (low  and  high  fidelity)  while 
demonstrating  several  numerical  and  physical  effects. 
The  numerical  consistency  of  the  solver  has  been 
demonstrated  by  conducting  in-depth  grid-sensitivity 


studies  ranging  from  0.25  mm  to  0.0625  mm,  or 
dnjdxi  ratios  of  0.588  to  2.352.  Successful  solutions 
were  obtained  with  resolution  levels  of  0.125  mm  for 
the  single-hole  spray  while  0.0625  mm  were  needed  for 
the  two-hole  spray.  The  additional  level  of  resolution 
needed  to  resolve  the  multi-hole  spray  arises  from 
spray’s  orientation  relative  to  the  numerical  grid. 
Because  the  grids  are  not  perfectly  aligned  in  this  case, 
numerical  dissipation  effects  may  emerge  causing 
unphysical  behavior  with  coarse  cell  sizes. 

Two  liquid  surrogates,  Surrogate-C  and  Modified 
Aachen  Surrogate,  were  used  to  model  JP-8  fuel.  Their 
selection  was  based  on  the  literature  for  gas-phase 
combustion  and  also  demonstrates  their  applicability  in 
the  liquid  phase.  Results  show  consistent  behavior 
among  the  sun'ogates.  This  was  demonstrated  in  the 
grid  sensitivity  analysis  as  both  surrogates  converged  to 
the  same  values  as  the  resolution  was  increased.  The 
behavior  of  the  surrogates  also  showed  good 
agreements  with  the  spray  measurements.  Moreover,  it 
was  established  that,  the  temporal  spray  breakup 
constant  required  adjustments  when  changing  from 
diesel  fuels  where  B^  —  7.  The  single-hole  axial  spray 
configuration  was  used  as  the  calibration  platform  for 
JP-8  fuel  resulting  in  calibration  of  the  spray  temporal 
breakup  constant:  —  10  (RANS),  and  B^  —  S  (LES). 
This  is  in  an  agreement  to  previously  developed 
physical  property  correlations  with  where  it  was 
shown  to  be  dependent  on  mass-density,  viscosity,  and 
surface  tension.  It  is  also  noted  that  when  using  high 
fidelity  models  the  RANS  constant  must  be  recalibrated 
with  experiments  as  it  was  reported  here. 

Low  fidelity  models  for  both  single  and  multi -hole 
simulations  show  good  comparisons  with  experiments 
in  terms  of  global  properties.  Because  of  the  coarseness 
of  RANS  approach  dissipation  will  dampen  out  the 
velocity  fluctuations  leading  to  inaccuracies  in  velocity 
magnitudes.  High  fidelity,  LES,  models  are  able  to 
capture  transients  and  are  a  useful  fool  to  provide 
insights  into  the  underlying  dynamics.  Using  the  LES 
model  the  interactions  between  two  adjacent  plumes, 
60°  apart,  were  found  to  be  weak  in  terms  of  affecting 
the  penetration  profiles.  This  was  concluded  after 
reporting  no  significant  change  in  liquid  penetration  for 
single  and  two-hole  nozzles  numerically  with  a  mean 
value  of  20  mm  for  both  cases  and  similar  vapor 
profiles.  The  experimentally  reported  decrease  in  liquid 
penetration  magnitude  from  single  to  two-hole  nozzles 
can  be  partly  due  to  non-linear  behavior  in  injected 
mass  quantities  (and  resulting  ROI)  that  was  observed 
in  Figure  2.  In  particular  for  asymmetrical  orifice 
configurations,  such  as  the  two-hole  60  degree  nozzle, 
local  Bernoulli  type  pressure  disturbances  near  the  sac 
volume  may  arise  and  result  in  non-uniform  flow 
distribution.  The  interaction  between  the  needle,  the 
fluid,  and  the  micro-channels  may  be  particularly 
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complex  in  this  case  resulting  in  the  observed  injection 
quantities  and  spray  behavior. 

Future  work  is  aimed  at  explaining  some  of  the 
reported  differences  by  conducting  a  deeper 
investigation  with  that  particular  nozzle  geometry. 
Detailed  nozzle  flow  calculations  can  be  carried  out  to 
extend  the  modeling  and  include  the  internal  flow 
dynamics,  fluid  structure  interaction  between  needle 
and  fuel,  and  the  emerging  spray.  In  terms  of  promoting 
enhanced  plume  to  plume  interactions  and  mixing, 
additional  work  will  also  address  the  effect  of  using 
cluster  nozzles.  Optimization  studies  with  respect  to  the 
number  of  cluster  used,  position,  and  operating 
conditions  will  be  carried  out  to  provide  cluster 
arrangements. 

In  summary,  some  of  the  major  findings  in  this 
study  are  as  follows: 

1.  Successful  grid  converged  solutions  can  be 
obtained  for  single-hole  and  multi-hole  sprays 
considering  the  challenges  imposed  by  the 
miss-alignment  between  the  spray  and  the  grid 
(for  multi-hole  spray).  The  methodology  for 
grid-converged  solution  was  established 
previously  [7]  and  has  been  applied  in  this 
paper  for  multi-hole  nozzles. 

2.  Similar  behavior  between  Surrogate-C  and  the 
modified  Aachen  Surrogate  was  established 
both  numerically  and  physically  with  good 
agreements  with  measurements.  Note  that  H/C 
ratio  is  maintained  fairly  uniform  across 
surrogate  mixtures  in  comparison  to  JP-8  with 
H/C  value  of  1.91. 

3.  Single-hole  LES  results  demonstrates  close 
agreements  with  measurements  with  JP-8 
sprays.  Liquid  and  vapor  penetration  length 
predictions  indicate  classical  high-pressure 
spray  features  beginning  with  a  highly 
transient  initial  stage  and  reaching  a  steady 
condition  for  liquid  penetration  only.  Vapor 
penetration  (and  mixing)  is  an  inherently 
dynamic  process. 

4.  Multihole  LES  results  show  the  effect  of  two 
adjacent  plumes  to  be  weak  at  60  degrees  with 
1mm  hole-to-hole  spacing  with  respect  to  the 
liquid  penetration  profiles.  The  decrease  in 
experimental  liquid  penetration  length  is  partly 
attributed  to  a  decrease  in  injected  mass  per 
hole  as  compared  to  the  single  hole,  and  to 
internal  flow  dynamics  and  interaction  with 
needle. 
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Figure  11.  Comparison  of  Single  Hole,  dynamic  structure  LES  model  (left 
column)  with  experiments  (right  column)  at  times:  33,  165,  495,  891,  1188  /rs. 
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Figure  12  1.  Evaporating  JP-8  spray,  single  hole  nozzle,  configuration  modeled  via  dynamic  structure  LES 
at  times:  33,  165,  495,  891,  1188  /rs.  Mixture  fraction  (1st  column).  Temperature  (2nd  column).  Turbulent 
Kinetic  Energy  (3rd  column),  Vorticity  (4th  column) 
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Figure  15.  Evaporating  JP-8  spray,  two  hole  nozzle,  configuration  modeled  via  dynamic  structure  LES  at 
times:  50,  100,  150,  250,  250  /rs.  Comparison  between  experiments  liquid  length  (left  column)  and 
instantaneous  LES  results  (right  columns)  -  Liquid  length  is  represented  by  particles  and  the  3D  vapor 

iso-surface  is  shown  on  grey. 
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